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Azimuthal correlations of high-pr charged hadron pairs from (di-)jet-fragmentation are studied 
at mid-rapidity in Au+Au collisions at y^s^=200 GeV. The distribution of jet-associated partner 
hadrons (1.0< pr <2.5 GeV/c) per trigger hadron (2.5< pr <4.0 GeV/c) is found to vary with 
collision centrality, both in shape and in yield, indicating a significant effect of the nuclear collision 
medium on the (di-)jet fragmentation process. 

PACS numbers; 25.75.Dw 



Energetic collisions between heavy ions at the Rela- This matter has been observed to strongly suppress the 
tivistic Heavy Ion Collider (RHIC) have been shown to yield of hadrons with large transverse momenta in cen- 
produce matter with extremely high energy density Jj. tral Au-|-Au collisions, compared to yields in p-l-p col- 
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lisions scaled by the number of binary nucleon-nucleon 
collisions 0, |M S ■ Such a suppression was predicted 
to result from energy loss of hard-scattered partons (light 
quarks and gluons) traversing the dense matter prior to 
forming the observed hadrons [EH- 1^ parton en- 
counters a sufficient amount of dense matter, the energy 
loss could strongly modify its fragmentation into jets of 
hadrons. 

Strong suppression of the away-side jet has been ob- 
served at RHIC Q. However, it is unclear at present 
how the lost energy is transported by the dense medium, 
and how the parton-medium interaction affects the frag- 
mentation process. Recently, there have been predictions 
that the coupling of jets to a strongly interacting medium 
may modifythe angular distribution and number of jet 
fragments EUSIIlIIlIIlIIlIIlIIi- Quarks from 
hard scattering processes may recombine with thermal 
quarks from the dense medium^, 0|. Co- moving radi- 
ated gluons may produce a "wake" in the medium, fur- 
ther increasing the number of quarks available for build- 
ing hadrons in the jet fragmentation process^30|. It 
has even been proposed that the energy deposited in the 
medium creates a shock wave around the propagating 
parton, thereby creatin g a "conical flow" akin to a sonic 
boom in a fluid [l^ . To investigate the transport 
of lost parton energy, the PHENIX experiment at RHIC 
measures azimuthally correlated hadrons arising from jet 
fragmentation as a function of centrality in Au-|-Au colli- 
sions. Such studies, in effect, use hard scattered partons 
as short- wavelength probes of the produced medium. 

The analysis presented in this Letter uses data from 
ys^=200 GeV Au-HAu collisions in the PHENIX Run-2 
data set. Charged particles are reconstructed in the cen- 
tral arms of PHENIX using drift chambers, each with az- 
imuthal coverage of 7r/2, and two layers of multi-wire pro- 
portional chambers with pad readout (PCI, PCS) |l7l| . 
Pattern recognition is based on a combinatorial Hough 
transform in the track bend plane, with the polar an- 
gle determined by PCI and the collision vertex along the 
beam direction |l8|. Particle momenta are measured with 
a resolution Sp/p = 0.7% © 1.0%p (GeV/c). To reject 
most background from albedo, conversions, and decays, 
a confirmation hit is required within a 2a matching win- 
dow in PCS The Au-|-Au event centrality is deter- 
mined using the PHENIX beam-beam counters (BBC) 
and zero-degree calorimeters (ZDC) T^. 

The traditional identification of jets through hadronic 
calorimetry and cluster algorithms is problematic in 
A-l-A collisions at RHIC, since low-energy jets (<10- 
20 GeV) are overwhelmed by other produced particles 
in the underlying event and high-energy jets are rela- 
tively rare at ys^=200 GeV. Instead, we study hard- 
scattered single partons and parton pairs through angu- 
lar correlations of high-pT hadron pairs. We examine 
the distribution of pairs over relative azimuthal angle, 
dN^-^ /d{A(j)), where A and B denote charged particles 



in the PHENIX pseudorapidity acceptance {\r]\ < 0.35) 
and in pr bins 2.5 GeV/c< p^ <4.0 GeV/c ("trigger"), 
and 1.0 GeV/c< <2.5 GeV/c ("partner"). Pairs 
from fragments of the same jet are expected to appear 
near A(f) ~ 0, while it indicates one hadron each 

from the outgoing hard-scattered parton pair (hereafter 
"dijet"). 

The PHENIX acceptance at central rapidity is non- 
uniform in azimuth. We correct for the shape of the 
acceptance in Acf) by constructing an area normalized 
correlation function, utilizing pairs from mixed events: 



C(A(/.) 



AB 



jYsiLiAcj^) rf(A0) 



(1) 



where 5^samc(A0) and Y^^^^{A(j)) are, respectively, the 
uncorrected yields of pairs in the same and in mixed 
events within each data sample. 
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FIG. 1: Correlation functions, C(A(j)), for pairs of charged 
hadrons with 2.5 GeV/c< p4 <4.0 GeV/c and 1.0 GeV/c< 
<2.5 GeV/c in different bins of collision centrality. The 
solid bands indicate the estimate of the background pair com- 
ponent (see text) within one unit of its systematic error. 

The correlation functions are shown in Fig. ^ folded 
into the range < A(f) < tt. For the most peripheral 
collisions (cf. Fig. ^), the correlation function shows 
two well-defined peaks centered at A(j) = and Ac/) — tt, 
which we can attribute to (di-)jet pairs. For more central 
collisions, there is a similar peak at A(j) = 0, a broader 
peak at A(j) — tt, and the apparent minima appear at 
A(f> < 7r/2. These features reflect a mixture of (di-)jet 
pairs and underlying events with particle flow along the 
reaction plane |20l |. 

In order to extract and examine the jet-induced pairs 
we analyze the pair distribution in the context of a 
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two-source model, assuming that each hadron can be 
attributed to either (i) a jet fragmentation source, or 
(ii) the underlying event. Additionally we assume that 
the (p distributions for A or B inclusive single particles, 
summed over both sources, have a shape proportional 
to (1 + 2(w^°"^) cos(2(0 - $Rp))) relative to azimuthal 
angle <i>Rp of the reaction plane of each event. All pairs 
which are not from the same jet or dijet fragmentation are 
termed background pairs, and are taken to have no angu- 
lar correlation beyond having their distributions respect 
the same reaction plane. In principle, contributions from 
resonance decays and global transverse momentum con- 
servation can also affect the distribution of background 
pairs, but we estimate these effects to be negligible for 
these pt ranges and the PHENIX rj acceptance. The 
distribution of background pairs over Acf) is then propor- 
tional to (1 -f 2(vfvf) cos(2A0)). 

Given a normalization, C(A0) can then be decom- 
posed into two pieces, one proportional to the distribu- 
tion of background pairs and another J( A0) proportional 
to that of the (di-)jet pairs: 

C(A0) = bo{l + 2{v^vi) cos(2A(/))) + J(A0) (2) 

We approximate (f^w^) — aud we measure 

(u^) and (f^) for each centrality and particle px bin 
through a standard reaction-plane analysis using the 
PHENIX BBC to reconstruct the reaction plane event- 
by-event. The large rapidity gap, A77 > 2.75, between the 
central arm acceptance and the BBC acceptance substan- 
tially reduces non-flow contributions to the measured V2 
values, particularly those arising from di-jets. The re- 
sults are shown in Tabled they are consistent with prior 
PHENIX V2 measurements j21i |. 

TABLE I: Anisotropy values for bins A (2.5 < pr < 4.0 
GeV/c) and B (1.0 < pt < 2.5 GeV/c), shown with statistical 
errors, and values of (j>Min (see text). The relative systematic 
errors on the anisotropics are estimated to be ±6% for the five 
most central samples and ±40% for the most peripheral sam- 
ple. The systematic errors on the V2's are dominated by the 
uncertainty in the correction for reaction plane resolution|2ll. 
and we assume them to be completely correlated between the 
two Pt bins in each centrality sample. 

Centrality (%) (vj) (vi} (t>Min (rad) 



0-5 


0.035 ± 0.001 


0.052 ± 0.007 


0.94 


5-10 


0.062 ± 0.001 


0.100 ± 0.005 


0.96 


10-20 


0.095 ± 0.0005 


0.144 ± 0.003 


0.98 


20-40 


0.146 ± 0.0004 


0.208 ± 0.003 


0.91 


40-60 


0.171 ± 0.001 


0.236 ± 0.006 


0.86 


60-90 


0.066 ± 0.001 


0.091 ± 0.004 


1.06 



The average level of the background, bo, can, in princi- 
ple, be fixed by making an assumption about the shape of 



the (di-)jet pair distribution. However, since we wish to 
measure the shape of the (di-)jet azimuthal correlations, 
we use a technique that requires no such a priori assump- 
tions. The simplest assumption allowing fcg to be fixed is 
that c?iVpf_-jjgj/(i(A(/)) is zero for at least one value of A^ 
(i.e. A(j)Min)- We refer to this as the ZYAM ("zero yield 
at minimum") assumption for the (di-)jet pair distribu- 
tion. The ZYAM condition is met by varying bo until the 
background component matches a functional fit to the 
correlation function at one point A(j)Mim as illustrated 
by the solid bands in Fig. ^ The systematic error on bo 
associated with this procedure (see Fig.|2J) was estimated 
by using a variety of functional forms that matched the 
data. 

A non-zero yield of (di-)jet pairs at A0Min would in- 
validate the ZYAM assumption and result in an overesti- 
mate of the value of 60. To verify that we are not making 
significant error in the normalization of the background, 
we have independently estimated the 60 values using the 
AB pair combinatorial rate, corrected for a slight bias in- 
troduced by mixing events of different multiplicity within 
the same centrality class These independent esti- 

mates are consistent with or even slightly higher than 
the ZYAM-determined values for bo, confirming that we 
are not significantly over-estimating the background lev- 
els. 

Once {V2V2) aud 60 are fixed, we can extract 
J(A(/)) and the fully corrected (di-)jet pairs distribution 
'^^(Df-)jct/^(^'?^)- construct the conditional yield 
distribution of jet-associated partners per trigger: 

1 dNf^t)j.t J(A0) N^^ 

d{A(j)) ^ C{A4>') d{A(t)') ^' 

Here, is the number of triggers and N^^ the total 
number of AB pairs in the event sample. Assuming that 
the pair efficiency is the product of the single particle effi- 
ciencies the trigger (A) efficiency cancels in Eq. |31 Thus, 
the ratio is corrected for acceptance and reconstruction 
efficiency |^ of the lower-pT B particles; the systematic 
error on this correction leads to a 10% uncertainty on the 
associated yields. 

The conditional yields of (di-)jet-induced partners per 
trigger are shown in Fig. [5] For the most peripheral event 
sample the (di-)jet associated yield distribution has an 
appearance we might expect from a normal (di-)jet frag- 
mentation process |23l |2J| : a well-defined near-side peak 
around Acj) = and a somewhat wider away-side peak 
around A<f> — tt. For more central event samples the 
shape of the near-side peak is essentially unchanged while 
the associated yield in the near-side peak increases, indi- 
cating some change in the fragmentation process . 

The much more dramatic change, however, is in the 
away-side peak which is considerably broader in all the 
event samples more central than 60%. In mid-central and 
central collisions there is a local minimum at A(j) = tt. 



5 




0.5 1 1.5 2 2.5 3 
A (j) (rad) 



0.5 1 1.5 2 2.5 3 
A (ji (rad) 



FIG. 2: Jet-pair distributions dN^_y-f^Jd{A4i) for diflfer- 
ent centralities, normalized per trigger particle. The shaded 
bands indicate the systematic error associated with the de- 
termination of A(j>Mi„. The dashed (solid) curves are the 
distributions that would result from increasing (decreasing) 
{V2V2) by one unit of the systematic error; the dotted curve 
would result from decreasing by two units. 



The existence of these local minima per se is not signif- 
icant once we take the systematic errors on (w^w^) into 
account (see below), but it is clear that the away-side 
peaks in all the more central samples have a very differ- 
ent shape than in the most peripheral sample. 

Given the dramatic results for the away-side peaks seen 
in Fig. 13 it is important to establish that they are not 
simply artifacts created by our method for background 
pair subtraction. If we relax the ZYAM assumption and 
lower bo slightly, the effect on any (di-)jet pair distribu- 
tion would essentially be to raise it by a constant, which 
would not change the presence of the local minima at 

Ac/) = TT. 

Changes to our estimate for (w^w^) can alter the shape 
of the (di-)jet distribution for some centrality samples, 
but the result of away-side broadening with centrality 
remains robust. The curves in Fig. [21 show the distribu- 
tions that would result if the (w^w^) products were arbi- 
trarily lowered by one and two units of their systematic 
error. With a two-unit shift the shape in the mid-central 
would no longer show significant local minima at A(j3 = tt. 
However, the widths of the away-side peaks are clearly 
still much greater than in the peripheral sample and the 
distributions in the two most central samples are hardly 



changed at all in shape. Even lower values of 
could be contemplated, but they would still not change 
the qualitative result of away-side broadening. And, such 



low {V2V2) values would also require a severe breakdown 
of the assumption (■y^f^) = indicating that 

these background pairs have a large, hitherto- unknown 
source of azimuthal anti-correlation. 

Convoluting the jet fragments' angles with respect to 
their parent par tons and the acoplanarity between the 
two partons (23 would yield a Gaussian-like shape in 
A0, possibly broadened through jet quenching (la l25l |. 
The observed shapes in the away-side peaks cannot result 
from such a convolution. 

We define the part of the A(j) distribution in 
\A(j)\ < A(/)Min as the "near-side" peak and \ A(j)\ > Acpuin 
as the "away-side" peak. Each peak is characterized by 
its yield of associated partners per trigger, and by its 
RMS width. We measure these for the full peak in the 
distribution over all values of Acj); the folded distributions 
over < A(p < tt shown here contain only half of each 
full peak's shape. These yields and widths are plotted in 
Fig. 121 for the different Au-|-Au centrality samples, along 
with the same quantities for 0-20% central d-|-Au colli- 
sions at GeV |2^. The yields and widths for 
the near- and away-side peaks in peripheral Au+Au col- 
lisions are consistent with those in d-|-Au collisions. The 
yields of both the near- and away-side peaks increase 
from peripheral to mid-central collisions, and then de- 
crease for the most central collisions. The near-side width 
is unchanged with centrality, while the away-side width 
increases substantially from the 60-90% sample to the 
40-60% sample and then remains constant with central- 
ity. 
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FIG. 3: (a) Associated yields for near- and away-side peaks in 
the jet pair distribution, and (b) widths (RMS) of the peaks in 
the full 0-2n distributions; plotted versus the mean number of 
participating nucleons for each event sample. Triangles show 
results from 0-20% central d-|-Au collisions at the same ys„„ 
[2^ . Bars show statistical errors, shaded bands systematic. 



In summary, we have presented correlations of high 
momentum charged hadron pairs as a function of col- 
lision centrality in Au-|-Au collisions. Utilizing a novel 
technique we extract the jet-induced hadron pair dis- 
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tributions and show that the dense medium formed in 
Au+Au colhsions at RHIC modifies jet fragmentation. 
In central and mid-central collisions the away-side an- 
gular distribution is significantly broadened relative to 
peripheral and d-|-Au collisions, and appears to be non- 
Gaussian. The shapes of the away-side A0 distributions 
for non-peripheral collisions are apparently not consis- 
tent with purely stochastic broadening of the peripheral 
Au-f Au away-side. However, the broadening and possi- 
ble changes in shape of the away-side jet are suggestive 
of recent theoretical predictions of dense medium effects 
on fragment distributions 0, 0, 0, |2^. The broad- 
ened shapes of the away-side distributions also imply that 
integration of the away-side peak in a narrow angular 
range around A0 = tt yields fewer associated partners 
in central collisions than in peripheral/d-|-Au collisions, 
as seen elsewhere [sl E^: but integrating over the entire 
broadened peak recovers the jet partners in the range 
1.0 GeV/c < pf < 2.5 GeV/c used here. Even though 
two-particle correlations do not allow for full reconstruc- 
tion of the jet fragmentation function, these data provide 
an entirely new way to probe the hot, dense medium 
formed in heavy ion collisions. 
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